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Abstract 

We investigate spherically symmetric non topological solitons in electro- 
dynamics with a scalar field self interaction [/ ~ IV'I taken from the complex 
signum-Gordon model. We find Q-balls for small absolute values of the total 
electric charge Q, and Q-shells when \Q\ is large enough. In both cases the 
charge density exactly vanishes outside certain compact region in the three 
dimensional space. The dependence of the total energy E of small Q-balls 
on the total electric charge has the form E ~ while in the case of 

very large Q-shells ^ ~ IQI'^/^. 



PACS: 11.27.-i-d,98.80.Cq, ll.lO.Lm 
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1 Introduction 



Several field-theoretic models predict existence of stationary finite energy excita- 
tions called Q-balls, see [1 , 2J for a review. Such objects arise as a highly nontriv- 
ial consequence of collective nonlinear dynamics of the fields. In contradistinction 
to topological defects which appear because of a frustration of the field interpo- 
lating between different degenerate vacua, the Q-balls exist in models with non 
degenerate vacuum. They must be taken into account in every thorough analysis 
of the dynamics of the fields. 

In a recent paper [[3]| we have shown that the Q-balls can be found also in the 
complex signum-Gordon model. They differ from other Q-balls by their strictly 
finite size: the vacuum field ip = outside a spherically symmetric Q-ball is 
reached at a finite radius R, and it is approached in a parabolic manner. Such 
parabolic behavior is the general feature of field theoretic models with V-shaped 
self interaction terms p]. It is a consequence of the fact that the field theoretic 
force density, which in the case of the complex signum-Gordon model with the 
complex scalar field ^ is given by —dU/dtp*, where U = Xl^pl and A > is the 
self coupling constant, does not vanish in the limit ^/^ — > 0. Therefore, for weak 
fields it is always stronger than the gradient force (~ Aijj) which alone would lead 
to an infinitely long tail in the spatial asymptotics of the field, and also stronger 
that the mass term (~ m^tjj) which together with the gradient force would lead to 
the well-known exponential asymptotics for the field. Heuristically, one can say 
that the parabolic behavior is due to the fact that the field resists with a finite force 
density even the slightest deviations from the vacuum configuration. The self 
interaction term f/dV'l) = Mi'l is V-shaped because its plot has the form of the 
inverted symmetric cone with the vacuum field ip = right at the tip. The models 
with V-shaped potentials are interesting for several reasons (4]. Let us mention 
here a scaling symmetry of the on-shell type, and the lack of a free field regime 
because the field equations can not be linearized around the vacuum solution. In 
case of smooth self interactions the second derivative of pertinent field potential 
at = gives a mass scale which is physically important in the vacuum sector. 
Defined in this manner mass scale is infinite if the potential is V-shaped. 

The motivation for the present paper comes from the idea that the finite vac- 
uum restoring force can perhaps win over the electrostatic repulsion for arbitrarily 
large electrically charged Q-balls. In order to address this question we extend 
the complex signum-Gordon model by including an Abelian gauge field min- 
imally coupled to the scalar field. Thus, the global U{1) symmetry is replaced by 
the local one. The considered model is a version of scalar electrodynamics. 
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It has turned out that in such scalar electrodynamics there exist electrically 
charged Q-balls provided that the absolute value of the total electric charge Q is 
small enough. For larger values of |(5| we find, rather surprisingly, Q-shells with 
completely empty interior. Ther is no upper bound on \Q\. Both Q-balls and Q- 
shells are static, spherically symmetric, and the charge density exactly vanishes 
outside certain finite, Q-dependent radius R. The scalar field approaches its vac- 
uum value ^/^ = in the parabolic manner, while the electric field of course has 
the standard Coulomb tail. The total energy is finite. Another very interesting 
finding is a zigzag in the plot of the total energy versus the total electric charge Q 
at intermediate values of Q, where the Q-balls transform into the Q-shells. 

In literature one can find several other versions of scalar electrodynamics with 
Q-balls, see [jSl [6l |7l [8]|. In all of them there exists a maximal value of the total 
electric charge the static Q-balls can have. Q-shells were not found. In our case 
the electric charge can be arbitrarily large, but there is the transition from Q-balls 
to Q-shells at certain intermediate value of the charge. 

The plan of our paper is as follows. In Section 2 we introduce the model 
and discuss field equations for the spherical solitons. The Q-ball solutions for 
small values of |Q| are presented in Section 3. Section 4 is devoted to the Q-shell 
solutions in the case of large values of \Q\. The transition from Q-balls to Q- 
shells at the intermediate values of IQI is studied in Section 5. Section 6 contains 
a summary and remarks. 



2 Preliminaries 

The Lagrangian of our model has the following form 

L = -\f,,F^'' + {D,i,yD^ij - (1) 

where if) is the complex scalar field in (3 + 1) -dimensional Minkowski space-time, 
A > is the self coupling constant, |^| is the modulus of ^, D^ij) = d^ip + iqA^ip, 
where g > is the electromagnetic coupling constant and is the U{1) gauge 
field. For the sake of convenience the fields tjj, A^, the space-time coordinates 
and the constants A, g are dimensionless. Of course, in physical applications they 
should be multiplied by certain dimensional constants. Euler-Lagrange equations 
give the Maxwell equations 

d,F^^ = f (2) 
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with the current density 

f = iq {ilj*d'''^ - d"-^*-^) - 2q^A''-^*-^, (3) 
and the scalar field equation 

By definition, the r.h.s. of Eq. (4) is equal to when = (see Section II in 
[f3l). The field potential f/d^M) = A|V^| can be regarded as the limit of the smooth 



potential Uai\^\) = \\J a + l^p when a -> 0+. 

We will consider the simplest static, spherically symmetric Q-balls with van- 
ishing magnetic field. The corresponding Ansatz has the form 

^ = e'<^^F{r), 1=0, Ao = Ao{r). (5) 

Here F(r) is a real function of the radial coordinate r and > is a constant real 
frequency. The spatial part of the Maxwell equations is trivially satisfied, while 
the Gauss law has the form 

AAo = 2q{u + qA^)F\ (6) 

where A is the three dimensional Laplacian. The scalar field equation is reduced 
to 

AF + (a; + gAo)2F=^signF, (7) 

where sign(O) = 0. Equations (6), (7) acquire a simpler form when we change a 
little bit the notation. With 

K = ^, G = V2qF, B = uj + 
v2 

we obtain the following equations 

Afi = fiG^ AG = -GB"^ + K signG. 

Note that the fields i?, G are gauge-invariant as opposed to Aq, if) and uj (the 
Ansatz (5) fixes the gauge only partially because the U{\) gauge transformations 
which depend only on time are still allowed). 
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Furthermore, the parameter k can be removed with the help of rescaUng 

B{r) ^ G{r) ^ K^/^G{K^'''r), r ^ k^^V, 

and therefore we may put k, = 1 without any loss of generality. After taking into 
account the spherical symmetry we finally obtain the following equations 

S" = --S' + SG'^ (8) 
r 

G" = - GB^ + signG, (9) 

r 

where ' denotes the derivative djdr. These equations are supplemented by two 
conditions. First, continuity of V^', V^o at r = implies the conditions 

S'(0) = = G'(0). (10) 

Second, because Eqs. (8), (9) together with the condition (10) are invariant under 
the transformations 

B{r) -B{r), G{r) -G{r), 
we may also assume that 

5(0) > 0, G(0) > 0. (11) 

We shall see that in fact 5(0) > 0. 

The total electric charge Q and the total energy E are given by the following 
formulas 

Q = [d^x f = -— Hdr r^BG^ (12) 
J q Jo 

{VAof + \DoiP\^ + Vr^iP + M " 



where 



E= drr^[B'' + G^ + B''G'' + 2K\G\]. (13) 
Jo 

As already signalled, we put k = 1. 

The Q-balls and Q-shells are represented by solutions of Eqs. (8), (9) which 
obey the conditions (10), (1 1) and have finite Q and E. Let us have a look at Eq. 
(8). It can be written in the integral form 



B'(r) = \ rdr'r'^B(r')G\r'), 
Jo 



Q = I drr^BG^ = (14) 



which implies that the function B{r) can not decrease when r increases because 
B'(0) > 0. Furthermore, we shall shortly see that in the case of Q-balls and Q- 
shells G{r) vanishes for all r > R, where R is the finite radius of the cloud of the 
electric charge. Therefore, B'{r) = Q/r'^ for r > R, where 

^ ^ qQ 

and in consequence, B(r) = [3 — Q/r for large r. The usual condition that the 
electrostatic potential Aq vanishes at the spatial infinity and the definition of the 
function B imply that (3 = uj. Hence, fox r > R 

Q 

B(r)=uo~=. (15) 
r 

It follows that Ao(r) = Q/{ATxr) for large r, as expected for the static spherically 
symmetric, spatially localized electric charge distribution. 

It is the well-known fact that Q-ball solutions in models with a global U{1) 
symmetry minimize the total energy E under the condition that the U(l) global 
charge is kept constant [[DUl. Also in the considered case of local f/(l) gauge 
symmetry our Q-balls and Q-shells minimize the total energy, but one has to as- 
sume not only that the total electric charge Q is constant but also that the Gauss 
law (8) holds. This latter condition follows from the fact that the Gauss law is in 
fact a constraint which is essential for the physical content of the electrodynam- 
ics. Therefore the variations 5B, 5G have to be interrelated in such a manner that 
the fields B + 6B, G + 5G also obey that constraint. It follows that in the linear 
approximation the admissible spherically symmetric variations obey the condition 

SB" = --6B' + G^6B + 2BG 6G. 
r 



Using this constraint one can show that 



fOO p 

5E = 2uj5Q -2 drr^ G" + 2G' /r + B^G - sign(G) 

— JO L 



5G. 



It is clear that the conditions 5E_ = 0, 5Q = Q imply Eq. (9). Moreover, it turns 
out that the second variation of the energy is positive because it can be written in 
the form 

a{6GY 



6 E = dr r 



{6B'Y + {6G'Y + {G6B + BdGy + 



[a + G2)3/2 



'Note that with our sign conventions the total electric charge Q is negative while Q > 0. We 
call Q the charge and Q the total electric charge. 
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(here we use the regularized version Ua{\i'\) of the field potential). Thus, the 
Q-ball and Q-shell solutions indeed give local minima of the energy. 

3 The Q-balls 

Qualitative analysis of solutions of the set of Eqs. (8), (9) can be accomplished 
with the help of the standard tool: a mechanical analogy. In the context of Q-balls 
such interpretation of field equations was used already in [9 |. Thus, we regard (8), 
(9) as Newton's equations of motion for a fictitious particle of unit mass moving in 
the (B, G')-plane. Because of the assumptions (11) we actually consider motions 
in the first quarter of the plane. The radial coordinate r plays the role of time. 
There are two forces acting on the particle: the friction force given by the two- 
vector {-2B'/r, -2G'/r) and the external force / = signG - B'^G). The 
external force does not have a potential, i.e., / 7^ — Vn with certain u. This fact 
reflects the different origins of Eq. (8) and (9): Eq. (8) is in fact the Gauss law for 
the electric field, hence it is a constraint, while Eq. (9) comes from the dynamical 
evolution equation for the scalar field. The B component of / (equal to BG^) is 
nonnegative. This means that the particle is constantly pushed in the B direction 
until (7 = 0. The G component of / vanishes along the hyperbola GB^ = 1 and 
along the straight line G = 0. It is negative above that hyperbola (a 'northern 
slope') and positive when < G < 1/B^ (a 'southern slope'). One may imagine 
that the particle moves in a curved valley with the bottom on the hyperbola, and 
that there is a '5-wind' which constantly pushes the particle in the positive B 
direction. 

The mechanical analogy presented above might seem relatively complicated, 
nevertheless it gives the right intuitions about the existence and the form of solu- 
tions of Eqs. (8), (9). For example, one can guess that perhaps the 'B-wind' can 
just push the particle along the bottom of the valley. This would mean that there 
exists a solution such that G{r)B'^{r) = 1. Indeed, it turns out that this is the case: 
the solution has the form 

B{r)=PoV^, G(r) = ^, 

r 

where /3q = 2/^/3, 70 = V3/2. Unfortunately, this solution is not interesting 
as far as the Q-balls are concerned because it is singular at r = and it has the 
infinite total energy and charge. 

Note that there are no Q-balls such that -B(O) = 0. In this case Eq. (8) and the 
condition -B'(O) = imply that B{r) = for all r > 0, i.e., that the particle moves 
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along the G-axis. Then, it follows from Eq. (9) that G{r) = G(0) + The 
electric charge density for this solution vanishes, but the total energy is infinite. 

The Q-ball solution corresponds to the trajectory of the particle which starts 
at the 'time' r = with the vanishing initial velocity -B'(O) = G'{0) = from a 
point (-B(O), G{Q)) lying on the 'northern slope', and reaches the S-axis at a finite 
'time' R, G{R) = 0. If at that 'time' also G'{R) = 0, then for r > i? the particle 
moves along the B-axis, asymptotically approaching the point {B, G) = {cu, 0), as 
it follows from formula (15). It is stopped by the friction force, while the external 
force / vanishes on the i?-axis. The condition G'{R) = means that the particle 
lands on the S-axis tangentially to it. This is possible because of the 'S-wind'. 
The trajectories of the type described above correspond to the simplest Q-balls. 
More complicated solutions can exist too, see the remark 4 in Section 6. 

We have found the solutions of Eqs. (8), (9) corresponding to the above dis- 
cussed trajectories numerically. For a given value of B(0) > we have adjusted 
the value of G{0) > 1/B^{0) until G{R) = = G'{R) for a certain value of R. 
Example of such a numerical solution is shown in Fig. 1. 




r 



Figure 1: Example of the Q-ball solution. For this solution B{0) =3.30, G{0) = 
0.49..., R =1.35..., B{R) =3.37..., Q = 0.05..., E =0.44..., a; = 3.41... . 



The numerical analysis is nicely supplemented by an analytical approximation 
which works in the case of very large values of (0) . It is based on the observation 
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that in this case the 'northern slope' is very steep and that the bottom of the valley 
lies very close to the i?-axis. Therefore, the particle very quickly reaches the 
bottom of the valley and climbs the 'southern slope'. The '5-wind' is relatively 
weak on that slope because it is proportional to G^, see the r.h.s. of Eq. (8). Thus, 
we may expect that the coordinate B of the particle is almost constant. 

Let us assume that B{r) is just constant, B{r) = B(0) for < r < R and that 
B{0) is very large. Then Eq. (8) is approximately satisfied: because GB^ ^ 1 
close to the bottom of the valley, the term BG^ is of the order B^^(0) and it 
asymptotically vanishes when B{0) oo. Equation (9) with constant B becomes 
ordinary linear differential equation. It can easily be solved by the standard meth- 
ods [fTOl in the intervals of r on which G has a constant sign. In the case G > 
the general solution which obeys the conditions G"(0) = has the form 

where A is a constant. The condition G{R) = gives 

R 

^^^52(0) sin(5(0)i?)' ^^^^ 

while the condition G'{R) = yields the following relation 

B{0)R = Co, (18) 

where cq is a solution of the equation 

tan Co = Cq. (19) 

Equation (19) has infinitely many solutions, but only for 

Co = 4.493... (20) 

the function G{r) has positive values in the whole interval < r < i? as required. 
Comparison with the numerical solutions shows that such approximate analytic 
solution is quite accurate already for B(0) = 4 and of course it becomes more 
accurate with increasing value of -B(O). 

The charge for the approximate solution (16) is given by the formula 
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(the integral over r is elementary). Formula (15) taken for r — R implies that 
CO = B{0) + Q/R, and with the help of (18), (21) 



.5.1/6 



Co Q 



-1/6 



+ 



Q 



5/6 



(22) 



The total energy E_ is calculated from formula (13). The integral is split into 
and J^. In the latter one only B' = Q/r^ does not vanish. Elementary 
integrations yield the following formula 



E 



5, 1/6 
6, 



12 



Q 



11/6 



(23) 



We see from formulas (18)-(23) that the radius R, the charge Q, and the total 
energy E_ decrease with increasing values of -B(O) - the Q-ball becomes smaller 
and smaller. Its charge density (~ BG^) exactly vanishes outside the radius R = 
(6cq/5)^/^ Q^^'^. Generally, the charge density has a maximum at the center (r = 
0), and it monotonically decreases towards zero, provided that Q is small enough. 
For larger values of Q the charge density is maximal at a finite radius r > 0. 

The approximate formulas (16) - (18) imply that the initial data (-B(O), G(0)) 
for the pertinent trajectories of the fictitious particle lie on the curve 



G{0)B\0) = 1 - 



cos Co 



5.603. 



(cos Co = —0.217...) which is plotted as the continuous line in Fig. 2. The thick 
dots denote the initial data determined numerically. It is clear that there is a very 
good agreement if B{0) > 4 what corresponds to Q < 0.018. 

The numerical investigations of the Q-balls for smaller values of -B(O) (that 
is S(0) < 4) have shown, rather unexpectedly, that the initial data (5(0), G(0)) 
move to smaller values of G(0), see Fig. 2. The largest value of G(0) is obtained 
for S(0) ^ 1.5. When B{0) is still decreased the points {B{0), G{0)) move to the 
'southern slope' of the valley. Close to the end of the G{0) > part of the dotted 
line in Fig. 2 the values of B{0) slightly increase. The end of the numerically 
determined Q-ball line has the coordinates B(0) = 1.317..., G{0) = (Q = 
9.753..., E = 73.237..., u = 4.546...). 

The charge Q and the total energy E_ of the Q-balls first grow as we move along 
the dotted line, but later they slightly decrease, see Section 5. The maximal values 
Q = 9.779... and E = 73.47... for the Q-balls are obtained for G(0) = 0.056..., 
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B{0) = 1.290.... On the other hand, the radius R always grows - its value at the 
end of the Q-ball line is R = 4.856.... The very interesting region of B{0) fa 1.3 
with small values of G(0) is discussed in detail in Section 5. 



G(0) 



0.25 



0.75 - 



0.5 



1 - 








2 



3 



4 



5 



6 



7 



B(0) 



Figure 2: The initial data for our numerical Q-ball and Q-shell solutions (repre- 
sented by the thick dots). The bottom of the valley is marked by the thin dotted 
line. The continuous curve is given by the equation G'(0)S^(0) = 5.603 obtained 
from the approximate analytic Q-ball solution (16). The thick dots on the line 
G{0) = correspond to the Q-shells, and those with G{0) > to the Q-balls. The 
thick dots represent the numerical solutions we have actually obtained - we do not 
suggest that the full set of initial data is discrete. 



The numerical solutions for the Q-balls show that for -B(O) ^ 1.3 and small 
values of G{0) the function G{r) and the charge density have maximal values at 
certain radius r > 0, while at r = there is a local minimum, see Fig. 3. In 
particular, for B{0) = 1.317..., G{0) = 0, that is at the end of the Q-ball line, the 
charge density exactly vanishes at r = 0. One may think that this Q-ball is empty 
at the center. Such interpretation immediately suggests the question whether there 
exist static Q-shells, for which the charge density does not vanish in a spherical 
shell with finite internal and external radiuses. It turns out that the answer is in the 
affirmative. 
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Figure 3: Example of a large Q-ball solution. The function G{r) has a local 
minimum at r = 0. For this solution B{0) = 1.2, G{0) = 0.221..., R = 4.76..., 
B{R) = 2.54..., Q = 9.52..., E = 71.2..., u; = 4.54... 



4 The Q-shells 

In the case of Q-shells the function G{r) exactly vanishes when r < tq and also 
when r > R> tq. Hence, the electric charge is localized inside the spherical shell 

given by the condition < r < R. The function B{r) has constant value equal 
to -B(ro) in the interval < r < tq, and the Coulomb form (15) when r > R. 
The corresponding trajectory of the particle from the mechanical analogy starts at 
certain 'time' r — ro > from the point (i?(ro), 0) on the S-axis, tangentially to 
the axis, moves across the valley and climbs the 'northern slope' until G' = at 
a 'time' ri > tq when it reaches a turning point. Next, it again crosses the valley 
and comes back to the _B-axis (tangentially to it) at the time R > ri. Solutions of 
Eqs. (8), (9) of that kind have been found numerically. An example is presented 
in Fig. 4. 

The initial data for the trajectory of the fictitious particle in the case of Q-shells 
have the coordinates (-B(O), 0), where B{0) = B{rQ). They lie on the horizontal 
axis in Fig. 2. The mechanical analogy suggests that the Q- shell solutions exist for 
arbitrarily large values of B{ro). The initial 'time' Tq (equal to the internal radius 
of the Q-shell) is determined by a balance between the friction term —2G'/r and 
the harmonic force —GB'^ in Eq. (9). The point is that the values of B grow 
monotonically with r because of the 'B-wind', and therefore the harmonic force 
is larger when the fictitious particle returns to the B-axis than during its flight to 
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the 'northern slope'. For this reason the returning particle would just cross the 
B-axis and move to the G < part of the {B, G')-plane were it not for the friction 
which can take away the right amount of energy from the particle (if tq is suitably 
adjusted), so that G' = exactly at the 'time' R when the particle reaches the 
B-axis. Thus, the friction term can not be neglected if the 'B-wind' is present. 

Note also that if we neglect the 'B-wind' (as in the case of large Q-shells 
discussed below) the friction would stop the particle before it returns to the B- 
axis. Therefore, as far as the Q-shell solutions are concerned the 'B-wind' and the 
friction go together - we may keep or neglect only both. 




Figure 4: Example of the Q-shell solution. The function B{r) is constant for 
r < To = 1, and almost constant in the layer 1 < r < 2 which Ues inside the 
Q-shell. For this solution tq = 1, B{0) = 1.452..., R = 5.06..., B{R) = 2.52..., 
Q = 10.27..., E - 78.00..., u; = 4.55.... 



In the case of large Q-shells significant simplifications appear and it is possible 
to find approximate analytic solutions. Because tq is large, the friction term in 
Eq. (9) may be neglected. Moreover, also B{ro) is assumed to be large. In this 
region of the {B, G)-plane the valley is narrow and its northern slope is very steep. 
Therefore, the fictitious particle will return to the B-axis rather quickly and the 
end point {B{R), 0) of its trajectory will be close to the starting point (-B(ro), 0). 
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Thus, we may assume that B{r) ^ -B(ro) in Eq. (9), i.e., that the 'B-wind' is 
not important too. Numerical solutions show that indeed B{r) is almost constant 
inside large Q-shells, see Fig. 5. 

B'(r) 
0.05 

0.04 

0.03 

0.02 

0.01 





100 100.3 100.6 100.9 101.2 101.5 



Figure 5: Example of the function B{r) with its derivative B'{r) for a large Q- 
shell. For this solution ro = 100, B{ro) = 3.71..., R = 101.69..., B{R) = 3.75..., 
B'{R) = 0.04..., Q = 501.8..., E = 7468.8..., u = 8.68.... At values of r larger 
than shown in the picture the function S(r) becomes constant and B'{r) decreases 
to zero. 




The simplified form of Eq. (9) is as follows (G > 0) 

G"{r) = -B^{ro)G{r) + 1. (24) 

Both the 'B-wind' and the friction have been neglected. Equation (24) has the 
following general solution 

G{r) = + Acos{B{ro)r + 6), (25) 

where A and S are constants. The Q-shell solution has to obey the following 
conditions 

G'(ro) = = G"(ro), G{R) = = G'{R). 
14 



In order to satisfy them we take 



^=-^^y S = -B{ro)ro, (26) 

and 

S(ro)(i?-ro) = 27r. (27) 

The relation (27) agrees very well with the numerical results already for B{ro) > 
3.0. For example, for -B(ro) = 3.083... we have obtained B{ro){R — ro) = 
6.261.... 

Relation (27) gives R provided we know ro and B(ro). B(ro) can be expressed 
by ro and Q. To this end we use the definition (14) in which B(r) = B{ro) and 
G{r) is given by formulas (25), (26). Elementary integration gives the following 
formula 

Q = 3nB-'^{ro){xl + 27rxo + ^n^/S - 5/2), (28) 

where 

Xo = ro-B(ro). 

It remains to determine ro. Equation (24) is not helpful here because it is 
invariant with respect to translations of r. We will determine ro by minimizing 
the total energy E_ under the condition that the charge Q_ has a fixed value. Here it 
is convenient to use another formula for the energy, namely 

E = QB{R) + % + r^{G'^ + 2\G\). (29) 

It follows from the definition (13) by applying Eq. (8) and integration by parts. 
Performing the (elementary) integration in formula (29) and eliminating R and 
B(ro) with the help of formulas (27), (28) we obtain a rather lengthy formula for 
E_ which contains only Q and xq. In the limit of large xq it can be written in the 
form 



1/6^5/6 



(30) 



where we have omitted terms of the order Xq or smaller. The minimum of E_ is 
obtained for 

^o = ^Q + Y- (31) 

Using this value of xq we find that in the limit of large Q_ 

^ ^ V3 2^/ V/^g^/^ (32) 
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i?(ro)=(^) Q'/'. (33) 

Then ro ~ Q^^^, and i? = ro + 27r/S(ro). Comparison with numerical results 
shows that formulas (30-33) are quite accurate already for Q = 150 - in the case 
of E_, -B(ro) they give numbers which differ from the corresponding numerical 
results by less than 1 percent, and for tq it is less than 7 percent. Of course the 
discrepancies decrease when we go to larger values of Q. 

5 The transition from Q-balls to Q-shells 

The transition from Q-balls to Q-shells occurs at = 9.753.... The largest Q- 
ball has the radius 7?* = 4.85..., and the initial data is G(0) = 0, B{0) = 1.317.... 
For such small values of B(0) the 'B-wind' can not be neglected. In consequence, 
the simple analytical approximations, successfully applied in the cases of small 
and large Q, do not work. We have studied the pertinent solutions of Eqs. (8), (9) 
numerically, concentrating mainly on the global characteristics: the energy E and 
the charge Q. 

Let us first briefly describe the general structure of the solutions. On the Q-ball 
side, the minimum of the function G at r = gradually deepens until (7(0) = 0. 
The particle from the mechanical analogy starts its motion from the 'southern 
slope'. The 'B-wind' carries it to the northern side of the valley (then G{r) in- 
creases), where the particle makes a turn - at that moment G{r) has the maximal 
value. Next, it again crosses the bottom of the valley and climbs the 'southern 
slope' until it reaches the G = axis at the 'time' r — R, tangentially to the 
axis (G'{R) — 0). For r > i? the particle moves along the B-axis. It slows 
down and approaches the point B{oo) = u in the Coulomb way given by for- 
mula (15). On the Q-shell side, the solutions have the shape presented in Fig. 4. 
Note the influence of the 'B-wind': B{R) = 2.52... is significantly larger than 
B{0) = B{ro) = 1.45.... 

Our main findings are presented in Fig. 6, where E and Q are plotted as the 
functions of the radius R (in the case of Q-shells it is the outer radius). It turns out 
that i? is a convenient independent variable because it monotonically grows when 
we change the initial data along the thick dotted line in Fig. 2 starting from the 
small Q-balls end (G(0) > 0, B{0) >> 1). Rather surprisingly, close to R^ the 
total energy and the charge of the large Q-balls decrease when R increases. The 
same is true for the small Q-shells with the outer radius slightly above i?*. Thus, 
the functions E(R), Q{R) have one local maximum and one local minimum. As 
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far as we can see from our numerical data the maxima occur at the same values 
of R. The same seems to be true for the minima. Note that the variations of 
K{R), Q{R) in this region of R are rather small. 

We knovv' that E_{.R) ~ ^'^^ for small Q-balls (i.e., when Q_ « QJ, and 
that E_{R) ~ Q^'^^ for large Q-shells (Q » QJ. Using the numerical results 
summarized in Fig. 6 we have constructed the plot of E_ versus Q for Q in the 
interval [9.73, 9.78]. It has the zigzag shape schematically shown in Fig. 7. 
The two spikes correspond to the maxima and minima of KiR) and Q{R). The 
presence of spikes in the plot of energy versus charge was found in other models 
with Q-balls long time ago, see, e.g., fTTl. 
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Figure 6: The plot of the energy E_ (the thick dots) and of the total charge Q_ (the 
triangles) versus the radius R in the transition region. In the case of Q-shells R is 
the outer radius. 



We expect that the spikes are infinitely narrow. Our argument is based on the 
formula 

dE _ dE/dR 
dQ ~ dQ/dR 
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If the maxima of E{R), Q{R) are located at Ri, R2, respectively, then close to 
them 



dE 
dR ~ 


>- 


R,f 


'd^E 
dR^ 


1 d'^E 
„ ^ 3 dR^ 


(R- 


Ri) + ... 


dQ _ 
dR ~ 




R2f 


'd^Q 
dR^ 


1 d^Q 
„ ^ 3 dR^ 

H2 


{R- 

R2 


R2) + ... 



where d^E_/dR^\R^, d^Q/dR^\R^ are different from zero (and <0). If Ri ^ R2 
the derivative dE_/dQ_ would become infinite at the point R = R2. The numerical 
results do not show any such increase of dE/ dQ, hence Ri = R2. 
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Figure 7: The schematic picture of the zigzag in the plot E_{Q)- The arrows 
point in the direction of growing R. The dot on the uppermost line corresponds 
to the solution which separates the Q-balls from the Q-shells (-B(O) = 1.317..., 
G(0) = 0). The upper spike corresponds to the maxima of E_ and Q, the lower 
one to the minima. 
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In consequence, dE/dQ is a continuous function of Rat Ri. This means that the 
direction tangent to the curve E(Q) changes infinitesimally when we infinitesi- 
mally cross the maximum. Completely analogous reasoning applies to the lower 
spike in Fig. 7 which corresponds to the minima of E_{R), QiR)- 

6 Summary and remarks 

1. Let us summarize our results. We have shown with the help of the mechanical 
analogy and of the numerical calculations that in the f/(l) gauged signum-Gordon 
model there exist compact, electrically charged non topological soli tons. They 
have the form of the Q-balls if the modulus of the total electric charge is smaller 
than Q^, = AnQ^/q, where = 9.753... and g > is the electromagnetic cou- 
pling constant. When |Q| > Q* we have found the Q-shells which are completely 
empty inside: the energy and the electric charge densities exactly vanish in a fi- 
nite ball around the center of the soliton. While electrically charged Q-balls have 
been found also in other models [S [6l |71 0|, the static Q-shells appear only in 
the gauged signum-Gordon model. Our feeling is that this is connected with the 
fact that the scalar field of the signum-Gordon model reaches its vacuum value 
= on a finite distance (in the parabolic manner), not only in the direction of 
the infinity but also in the direction of the origin. Of course, also the electrostatic 
repulsion is an important ingredient - it prevents the Q-shell form a collapse. 

Furthermore, we have found the approximate analytic solutions of the field 
equations which agree very well with the numerical results in the limits of very 
small and very large absolute values of the electric charge. This is especially 
important in the latter case because the purely numerical approach is not sufficient 
in order to convincigly tackle the limit \Q\ — > oo. 

2. It is interesting to see in detail the change of the exponent from 5/6 to 
7/6 in the dependence of the energy E_ on the charge Q when passing from the 
asymptotic regions of very small to very large charges. Note that the Figure 7 
reminds a first order phase transition in, e.g., van der Waals gas, lfT2ll . We have 
not found a good analytic approximation for solutions in the intermediate region. 

3. The formula given at the end of Section 2 shows that the second variation 
of the energy 6'^E_ is positive. Therefore the Q-balls and the Q-shells are stable 
against small radial perturbations. In general however the question of their stabil- 
ity is rather complex one if one considers also nonsymmetric perturbations. First, 
they can perhaps decay by emitting small Q-balls and/or a scalar and electromag- 
netic radiation. We have checked on our numerical solutions that a split into two 
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objects of the charge (5/2 is energetically forbidden only for small Q-balls {Q < 
1 .62...). This of course does not prove that such a decay of large Q-balls in a finite 
time is dynamically possible. Furthermore, the evolution of the perturbed soli- 
tons can be very nontrivial because of complicated electromagnetic interactions 
of currents and charges. Note also that the usual criterion for the global stability 
consisting in excluding the possibility of evaporation by emission of quanta of the 
scalar field ^ is hard to apply in our case because the signum-Gordon model can 
not be linearized if the scalar field is close to the vacuum field. Therefore, an esti- 
mate of the rest mass of the quanta is not available - it is a rather difficult problem 
which belongs to non perturbative quantum field theory. 

4. We have investigated only the simplest (elementary) Q-balls and Q-shells. 
The mechanical interpretation of Eqs. (8), (9) suggests that there also exist solu- 
tions such that the charge density has several maxima. In this case the classical 
particle travels along the valley for a longer time climbing the 'southern slope' 
more than once until it finally settles on the B-axis. Also numerical investigations 
have shown solutions of this type. The relation E_{Q) for such complex Q-balls 
and Q-shells is not known as yet. Another open problem is their stability against 
a decay into the elementary Q-balls and Q-shells. 
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